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Abstract: We investigate the Dyson-Schwinger equations for the gluon and ghost propagators 
and the ghost-gluon vertex of Landau-gauge gluodynamics in two dimensions. While this simplifies 
some aspects of the calculations as compared to three and four dimensions, new complications arise 
due to a mixing of different momentum regimes. As a result, the solutions for the propagators are 
more sensitive to changes in the three-point functions and the ansatze used for them at the leading 
order in a vertex expansion. Here, we therefore go beyond this common truncation by including the 
ghost-gluon vertex self-consistently for the first time, while using models for the three-gluon vertex 
which reproduce the known infrared asymptotics and the zeros at intermediate momenta as observed 
on the lattice. A separate computation of the three-gluon vertex from the results is used to confirm 
the stability of this behavior a posteriori. We also present further arguments for the absence of the 
decoupling solution in two dimensions. Finally, we show how in general the infrared exponent k of the 
scaling solutions in two, three and four dimensions can be changed by allowing an angle dependence 
and thus an essential singularity of the ghost-gluon vertex in the infrared. 
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1 Introduction 



Quantum chromodynamics (QCD) is the theory of quarks and gluons. It is well understood in the 
perturbative regime, but phenomena like confinement and chiral symmetry breaking are intrinsically 
non-perturbative and their investigation therefore requires adequate methods. Some of these effects 
are believed to be present already for the gluonic sector alone, i. e., Yang-Mills theory, on which we will 
focus here. One approach to understand the non-perturbative properties is based on the correlation 
functions, which are also used as input in many phenomcnological applications. However, correlation 
functions are in general gauge-dependent. 

A common gauge choice is the Landau gauge, for which propagators and vertices have been cal- 
culated with several methods, see, for instance, [1-25], summarized in the recent review [26]. For 
various reasons investigations were extended from four also to three [3, 4, 18, 27-30] and two dimen- 
sions [3, 4, 28, 31-34]. One of the main motivations for this is that lattice calculations become much 
cheaper in lower dimensions so one can more easily reach low momenta. An approach complementary 
to lattice calculations is based on functional equations [35] like functional renormalization group equa- 
tions, see, e. g., [36-39], or Dyson-Schwinger equations (DSEs), see, e. g., [26, 40-45], for which large 
scale separations are less problematic. However, since they consist of an infinite tower of equations, 
truncations are required. 

The situation in four dimensions is interesting, because with continuum methods one finds two 
types of solutions, while lattice methods yield only one. The two types are called decoupling and 
scaling solutions. The former [11, 13, 16, 17, 20, 21] is characterized by a gluon propagator that 
becomes finite at zero momentum while the ghost propagator behaves like a massless particle. On the 
other hand, the scaling solution [1-4, 20, 21] features an infrared (IR) vanishing gluon propagator and 
an IR enhanced ghost propagator. On the lattice only the decoupling type of solution is seen in four 
and three dimensions. Currently there does not exist a consensus in the community if only one solution 
is physical or both types are valid and possibly correspond to different non-perturbative completions 
of the Landau gauge, see [26] for a compilation of the current state of affairs. The situation in three 
dimensions is essentially the same. 

On the other hand, the two dimensional case is different: Lattice calculations [31, 33] seem not to 
find the decoupling type solution, but the results resemble more the scaling type solution. However, 
studies in the strong coupling limit show that the situation is not completely clear, since the analytically 
known scaling relation for the IR exponents is not fulfilled for j3 = and the coupling does not become 
a constant [46] . The results for the dressing functions at finite f3 fit reasonably well with results from 
analytic analyses using functional methods [3, 4, 28, 34], but an agreement between the IR exponents 
extracted from the lattice data [31, 33, 46, 47] is not reached. Results obtained within the Gribov- 
Zwanziger framework [32] corroborate the non-existence of the decoupling type solution as does Ref. 
[34]. In the latter also a bound on the gluon propagator at zero momentum incompatible with this 
type of solution was derived. 

In the present work we want to extend the currently available analytic results from DSEs in two 
dimensions [3, 4, 28, 34] and present the first solutions at all momenta. We confirm the absence of the 
decoupling solution from Dyson-Schwinger equations in agreement with lattice and prior results also 
when taking the full momentum range into account. Although the underlying equations seem easier 
to solve in two dimensions than in four dimensions, for example, because no renormalization has to be 
performed, we find that they have their own intricacies. The reason is that in two dimensions different 
momentum regions influence each other. Especially the mid-momentum behavior, which is where the 
truncations of DSEs are most influential, is very important in order to obtain correct results. This can 
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be uniquely traced back to the subtle cancellations required for the scaling solution in two dimensions. 
As a consequence, it is non-trivial in two dimensions to find a viable truncation which works over the 
whole momentum range for the deep infrared to the ultraviolet. As a result we will discuss several 
model ansatze for the vertices. As a final step we also include the ghost-gluon vertex to the set of 
equations to be solved numerically. In this system of equations the only undetermined quantity is the 
three-gluon vertex, which we adjust such as to obtain the correct ghost UV behavior. 

We will explain the setup of the equations and fix our notation in Sec. 2. Then we will consider 
the ghost DSE numerically in Sec. 3. In Sec. 4 we present analytic results and the coupled system 
of the propagators DSEs is investigated numerically in Sec. 5. In Sec. 6 we present the results of 
the system of propagators and ghost-gluon vertex. We conclude with a summary in Sec. 7. Two 
appendices contain details on the integral kernels and results for the three-gluon vertex. 



2 Dyson-Schwinger equations of two-dimensional Yang-Mills theory 

The Euclidean Langrangian of Yang-Mills theory fixed to the Landau gauge reads 

where F^ v is a component of the field strength tensor = F® v T a with the Hermitian generators T a 
of the gauge group and M ab is the Faddeev-Popov operator: 

F£ v := d^Al - d v Al +igf ahc AlAl, (2.2) 

M ab ._ _ 5 abQ2 _ gfabeg^ (3 3) 

The DSEs can be derived from this Lagrangian with standard methods, see, for example, [40, 42, 44]. 
Since vertex DSEs are already rather complex the program DoFun [44, 48] was used for their derivation. 
The propagators of the ghost and gluon fields are given by 

D gh (p 2 ) := -^jp, D gl (p 2 ) := P^p)^ 1 , (2-4) 

respectively, where the transverse projector is P^ v (p) — 5^ v — p^p^/p 2 . Their full DSEs are given in 
Fig. 1. In the IR we will use the parametrization 

G IR {p 2 )^B-{p 2 ) 5 o\ Z IR (p 2 ) = A ■ (p 2 ) Sa ' (2.5) 

for the dressing functions, where S g h and 5 g i are the so-called IR exponents which describe the IR 
behavior of the dressings qualitatively. For the scaling solution they are related by [3, 4] 

25 gh + 5 gl =^^. (2.6) 

Consequently, in the case of the scaling solution the IR behavior of both dressing functions can be 
described by one variable only which we will denote by k := —5 g h- 

For the two-point DSEs we adopt a standard truncation scheme from four dimensions where all 
diagrams involving the bare four-gluon vertex are dropped. Such a truncation conserves the leading 
IR and UV behavior. In order to obtain a scalar equation from the gluon DSE we project it with the 
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Figure 1. Full two-point Dyson-Schwinger equations of Landau gauge Yang-Mills theory. All internal propa- 
gators are dressed. Thick blobs denote dressed vertices. Wiggly lines are gluons, dashed ones ghosts. 



transverse projector, yielding 
1 



G(p 2 ) 



l + N c g 2 Z(q')G{{ P + qY)K G {p, q )T Acc {q-p + q lP ) (2.7) 



' =1 + N c g 2 I G(q 2 )G((p + q) 2 )K 9 z h (p, q)T A5c (p;p + q, q) 



Z(P 2 ) 

+ N c g 2 ( Z{q 2 )Z{{p + q) 2 )K^(p,q)T A \p,q,p + q). (2.8) 



The quantities T Acc and T A are dressing functions of the ghost-gluon and three-gluon vertices, re- 
spectively, and J stands for J d 2 q/(2i{) 2 . The kernels Kq, K% h and are given explicitly in 
Appendix A. 

The dressed ghost-gluon vertex is described by two dressing functions 

Tf c > abc (k;p,q) :=igr bc {P^{k)p v Df c {k ] p,q) + k^Df cc (k-p,q)). (2.9) 

The basis tensors have been chosen such that D Acc {k\p 1 q) and D Acc (k;p, q) are the purely transverse 
and longitudinal dressing functions, respectively. After transverse projection then only D Acc (k;p,q) 
contributes and T Acc (k;p, q) from eqs. (2.7) and (2.8) can be identified with Df cc (k;p,q). 1 The bare 
vertex is ig f ahc Ppi- 

The three-gluon vertex has in general 14 Lorentz tensors. Only four of them are transverse and 
thus contributing in the Landau gauge. Furthermore we only consider the color antisymmetric part. 
As part of our truncation we only take into account the tree-level tensor: 

r^ p abc b, Q,r)=ig f abc (g MV (q ~ p) p + g vp (r - q) M + g pfl (p - r)„) D A * (p, q, r). (2.10) 

Thus T A3 (p, q, r) in eqs. (2.7) and (2.8) is identical to D A3 (p,q,r). 

Specific expressions for the dressing functions of the vertices will be given below. 

the notation of [4] we have D^ cc (k; p, q) = A(k; p, q) and D^ cc (k; p,q) = B(k; p, q) + A(k; p, q) p ■ k/k 2 . 



- 4 - 



3 The ghost Dyson-Schwinger equation 

This DSE is the simplest one because it has only one integral with a relative simple structure (compared 
to the gluon loop of the gluon DSE) and is automatically UV finite. Consequently it is in general 
relatively easy to study numerically using input for the gluon propagator and the ghost-gluon vertex, 
see, for example, [16]. In the two-dimensional case the ghost DSE provides additional important 
information. First we investigate the existence of various types of solutions, then the influence of the 
mid-momentum regime on the UV behavior of the ghost propagator is scrutinized. In this section we 
use a bare ghost-gluon vertex and various ansatze for the gluon dressing function. 

3.1 Existence of decoupling and scaling solutions 

In functional equations the choice between decoupling or scaling solutions is realized via the ghost 
Dyson-Schwinger equation. Employing the subtracted DSE, 



we have to specify the value of the ghost dressing function at pq. Most conveniently one may choose 
Po = 0. For a finite value of 1/G(0) a decoupling solution is obtained, whereas a scaling solution 
corresponds to 1/G(0) = 0. When solving the coupled system of the two-point DSEs the gluon 
propagator then automatically becomes of the decoupling or scaling type in four dimensions [2, 4, 21]. 

The unsubtracted ghost DSE in two dimensions was considered in Ref. [34], where it was found 
that in order to avoid IR divergences at p 2 = the gluon propagator has to vanish at zero momentum. 
Here we will corroborate these results numerically. Therefor we use various ansatze for the gluon 
dressing function, which are depicted in Fig. 2: 



For now the parameter c g i is set to one. The IRE exponent S g i is either 1 for the decoupling type 
ansatz or 1.4 for the scaling type ansatz. 

The subtracted ghost DSE in Eq. (3.1) was solved using a standard fixed point iteration. The 
resulting ghost propagator dressing functions are shown in figs. 3 and 4. For the scaling type ansatz 
the expected power law is perfectly obeyed in the IR, while for the decoupling type ansatz no valid 
solution is found. First of all the obtained ghost dressing function, depicted in Fig. 4, does not become 
constant in the IR but diverges. Fitting the IR exponent locally, i. e., extracting the IR exponent from 
two adjacent points, shows that the exponent does not settle to a constant value but keeps decreasing 
slowly. Secondly the ghost dressing depends on the value of the used IR cutoff, as shown in Fig. 4, 
hinting at the occurrence of IR divergences. Trying several variations of the gluon dressing function 
ansatz could not remove this problem as expected from the analytic considerations of Ref. [34]. For 
the scaling type solution the results for the same IR cutoffs as used for the decoupling type solution 
lie on top of each other. 

Note that both the analytic [34] and the numeric arguments are based on the assumption that the 
ghost-gluon vertex is constant in the IR and does not deviate much from the tree-level at intermediate 
momenta. From lattice calculations alone, however, we can strictly speaking not infer that for a 
decoupling solution, because what we see there seems to be a scaling solution. Also the often cited 
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Figure 2. Input gluon propagator dressing functions Figure 3. Ghost propagator dressing function re- 
of decoupling and scaling types from Eq. (3.2) with suiting from the scaling type gluon propagator. 
c g i — 1 and k = 1 (decoupling) and n = 1.4 (scal- 
ing). This plot and the following were created with 
Mathematica [49] 

G(p 2 ) 

0.75 
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Figure 4. Ghost propagator dressing functions re- 
sulting from the decoupling type gluon propagator. 
The various curves correspond to different IR cutoffs. 

Taylor argument [50] does not forbid a ghost-gluon vertex that is IR vanishing. Of course such a 
behavior is unlikely and unexpected, but one can construct a ghost-gluon vertex that allows a solution 
for the ghost equation also in this case. We also tried the opposite way, namely to solve the ghost- 
gluon vertex DSE using decoupling type ansatze for the propagators. However, in this case we could 
not find a solution. This further corroborates the non-existence of a decoupling type solution in two 
dimensions. 

3.2 Influence of the mid-momentum regime on the ghost's UV behavior 

Another important lesson from the ghost DSE alone concerns the influence of different momentum 
regions on each other. In four dimensions such an influence is subleading, mainly due to renormaliza- 
tion. In three dimensions, this is a weak quantitative effect [27]. For example, the UV behavior can be 
determined self-consistently. In two dimensions this is not the case. This can already be inferred from 
a purely perturbative investigation based on bare propagators which fails because of IR divergences 
which can only be remedied by taking into account a non-trivial behavior at momenta below the UV 
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Figure 5. Input gluon propagator dressing functions 
given by Eq. (3.2) and c gl = 1, 0.7, 0.4 (solid/red, 
dashed/green, dotted/blue). 



Figure 6. Ghost propagator dressing functions re- 
sulting from the gluon propagator ansatze given in 
Fig. 5 and a bare ghost-gluon vertex. Solid/red, 
dashed/green and dotted/blue lines correspond to 
c gl = 1, 0.7, 0.4. 



regime. We explicitly demonstrate the influence of the mid-momentum regime by studying the ghost 
DSE with the gluon dressing function given in Eq. (3.2) where we vary the parameter c g i which modi- 
fies the height of the bump in the gluon dressing function. The ansatze for c g i = 0.4, 0.7, 1 are shown 
in Fig. 5 and the resulting ghost dressing functions in Fig. 6. It is clearly visible that the height of 
the bump is correlated with the value of the ghost dressing function in the UV: The higher the bump 
the lower the ghost dressing becomes in the UV. This observation will be important later in Sec. 5 
when we solve the DSEs of both propagators simultaneously. Note that the IR is not affected by the 
mid-momentum behavior. 



4 The IR solutions revisited: Analytic considerations 
4.1 The solution k = 

For the scaling type solution one can determine the qualitative behavior in the deep IR by calculating 
the IR exponents which can be done analytically [1, 3, 4]. Therefor one calculates the IR leading 
diagrams in the two-point DSEs, viz. the ghost loop in the gluon DSE and the loop of the ghost 
DSE, see Fig. 1, using the power law ansatze for the propagator dressing functions given in Eq. (2.5) 
together with the scaling relation Eq. (2.6): 

Z IR (p 2 ) = A-(p 2 ) 2K+ ^/ 2 , G IR {p 2 ) = B-(p 2 )- K . (4.1) 

For the ghost-gluon vertex usually a bare vertex is employed but the same results for k hold for all 
ghost-gluon vertices with a regular IR limit [4]. The IR value of the coupling, however, depends on 
the IR value of the vertex. The calculation of k in d dimensions boils down to solving the following 
equation [3, 4]: 

sin(7r n)T(d/2 - «)r(«)r(l + d/2 + k) 
2(d - 1) sin(7r(d/2 - 2 K ))T(d - 2«;)r(2K)r(l + k) ' ^ ' ' 

For two dimensions two solutions have been cited in Ref. [3]: K\ — 0.2, which is close to values 
extracted from lattice data [31, 46, 47], and K2 — 0. Note that the IR exponents of the propagators for 
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the second solution are S g h = and S g i = 1 which are exactly the IR exponents of the decoupling type 
solutions. Hence the two-dimensional case is special in that a possible decoupling solution would also 
respect the scaling relation Eq. (2.6). We will now present further evidence against the decoupling 
type solution with K2 in addition to those presented in [34]. 

Let us scrutinize Eq. (4.2) for d = 2 and k = 0. In Ref. [3] it was argued that for d = 2 + e a 
solution k = + e exists. This is illustrated in Fig. 7, where the left-hand side of Eq. (2.6) is plotted 
for d = 1.9, 2, 2.1. As d approaches 2 from either side, there are intersection points with unity which 
indeed approach n = 0, such that there seems to be a solution also for d = 2. On the other hand, we 
can take the limit d — > 2 on the left-hand side analytically: 

1 + " 1. (4.3) 



2 -4k 



The only solution to this is k = 0.2. The reason for this discrepancy is that d = 2 and k = 
corresponds to a multi-valued point and a wide range of possible values can be obtained. For example, 
we could approach the point by taking the limit e — > and d = 2 + 2e and k = e. Then the left-hand 
side of Eq. (4.2) is one and k = is a solution. Note that in four dimensions the same arguments 
apply to the solution k = 1. In three dimensions the mathematical existence of two solutions is 
unequivocal. Varying the number of dimensions continuously one sees that there are two branches of 
solutions [3, 27]. In four dimensions only one of those two solutions is found to be realized in numeric 
calculations [51], but in three dimensions solutions for both cases have been found [27]. Nonetheless, it 
may be hypothesized that only one branch is physical, namely the one on which n = 0.595353 in four 
dimensions and k = 0.2 in two dimensions lie. The presence of the second solution in three dimensions 
could be a truncation artifact, but this requires further investigations. 

In summary, the infrared exponent k 2 requires an additional prescription for how to define a 
two-dimensional solution in a particular limit d — > 2 suitably combined with n(d) — ¥ 0. This makes 
its existence scheme dependent and it seems at least questionable whether such a solution is realized. 
Since this solution coincides with the decoupling solution, such a type of solution seems, from this 
point of view, unlikely to exist in two dimensions. 

4.2 Scaling solutions with k ^ 0.2 

Although lattice calculations seem to obtain a scaling type solution, it is not yet settled what the 
numerical value of k is. First results corresponded to values of about 0.15 [31], but recent calculations 
on larger lattices seem to favor a value of 0.225 [47]. From functional equations, however, the value 
k = 0.2 is well established [3, 4]. In the calculation of n only the IR parts of the propagators and the 
ghost-gluon vertex enter. Hence modifications in the mid-momentum and UV regime are irrelevant as 
are any details of the three-gluon vertex as long as it is not more IR divergent than (p 2 )~ 3fi +( d_4 )/ 2 
[9, 28]. As far as the propagators are concerned only the combination AB 2 of the coefficients of the 
IR power laws from Eq. (2.5) enters and is calculated together with k. The only quantity which has to 
be set by hand is the ghost-gluon vertex in the IR. A class of regular vertex dressings was investigated 
in Ref. [4], where it was shown that they all result in the same k. Regularity means here that the 
ghost-gluon vertex dressing function has a unique IR limit when all momenta go to zero. Thus the 
only way to obtain a different value for k is to allow non-regular vertex dressings, i. e., the IR limit 
depends on the angle between two momenta. 

To implement such an angle dependence we employ the ansatz 

D^ d {k-p,q) = l + a{co^) 2 (4.4) 
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Figure 7. Left- and right hand sides of Eq. (4.2). The dashed, solid and dotted lines corresponds to d = 
1.9, 2, 2.1, respectively. As d approaches 2, a solution seems to be k = 0, but actually the point d — 2, k — 
is multi-valued. The blob at k = 0.2 represents the standard solution. 



for the transverse dressing function of the ghost-gluon vertex, where ip is the angle between the two 
ghost momenta: 

7 2 2 2 

/ \ P ' Q k — p — q 

COS W = ] -rr- = — rr-j • 4.5) 

mm 2 \p\m 

The angle is not well defined in the zero-momentum limit. 

The ansatz (4.4) is motivated by the ghost-antighost symmetry in Landau gauge and corresponds 
to the first two terms in a Fourier series of more generally possible angle dependences [42] . The values 
for k and the IR fixed-point value for the coupling, given by [2, 4, 28, 52] 

a{0):=^-AB\ (4.6) 

can then be calculated using this ansatz. The results are depicted in Fig. 8. Figs. 9 and 10 show 
results for three and four dimensions, respectively, where the same arguments as in two dimensions 
apply. 

Finally we want to make a comment on the dependence of k on the employed projection of the 
gluon DSE. In principle any projection should yield the same value for k, but since we have to use 
a truncated system of equations it actually does depend on the projection operator. For example, 
the Brown-Pennington projector leads to k = 0.224745. However, any non-transverse projection leads 
to contributions from the longitudinal part of the ghost-gluon vertex. If we had the exact ghost- 
gluon vertex, this would be fine and we would get the same value of k for all projections. Since we 
do not have the exact vertex, we will employ the transverse projection only. This disposes of the 
following ambiguity as shown in Ref. [4]: In general there exists a second gauge parameter which 
interpolates between the normal Faddeev-Popov action of linear covariant gauges and the ghost anti- 
ghost symmetric gauges. However, in the Landau gauge this parameter is irrelevant and thus any 
quantity depending on it is ambiguous. A transverse projector gets rid of all such terms. See also Ref. 
[53] for a short discussion of this issue. 



or(0) 

0.4 r 



Figure 8. Left: n for different values of a. Right: The corresponding values of a(0). Both plots are for two 
dimensions. 



or(0) 
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Figure 9. £e/f: k for different values of a. Right: The corresponding values of a(0). Both plots are for three 
dimensions and only one solution branch is shown. 



5 The coupled system of two-point DSEs 

As a next step we solve the DSEs for the ghost and gluon two-point functions simultaneously. This 
system contains as a new quantity the three-gluon vertex for which we need an ansatz. An expression 
similar to the one used in four dimensions [51, 54] turns out to be insufficient, because the gluon loop 
will then dominate the ghost loop at intermediate momenta and the gluon dressing function turns 
negative. This was already observed in three dimensions [27] and motivated another construction for 
the three-gluon vertex. A similar one for two dimensions is 

V ans (p,q,r) - ^ Z (p2)Z(q2)Z(k2) J' {bA) 

This expression becomes 1 for large momenta and is suppressed for low ones. The exponent a controls 
the strength of the suppression. In order to find a solution it must not be too small. However, it 
is known that the three-gluon vertex is not IR suppressed but actually IR enhanced in two [28, 31], 
three [18, 28] and four dimensions [9, 18]. In two dimensions, an approximate agreement between 
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Figure 10. Left: k for different values of a. Right: The corresponding values of a(0). For a > 1.3 the coupling 
becomes negative and is not shown. Both plots are for four dimensions. 



the prediction from functional methods [28] and the lattice results [31] is found, while in higher 
dimensions the lattice results [18] do not yet penetrate very far into the infrared. Furthermore, lattice 
calculations showed that the three-gluon vertex changes its sign at intermediate momenta [18, 31]. 
For two dimensions we show in Appendix B that the DSE of the three-gluon vertex reproduces this 
zero crossing. Based on this behavior we adopt the following ansatz for the dressing function of the 
three-gluon vertex: 



+ h IR (p 2 + q 2 + r 2 Y 



(5.2) 



p 2 + q 2 + r 2 + h IM h- 2 

The first term ensures the correct UV behavior and the second term implements the IR divergence 
[28] . Note that the parameter hjn has to be negative in order to reproduce the zero crossing. h]R can 
be obtained by a fit to lattice data. However, in order to obtain a solution we have to choose a much 
higher absolute value in order to compensate for dropping the two-loop diagrams in our truncation. 

As in four dimensions the gluon DSE contains spurious divergences. They are the remnants of 
the quadratic divergences in four dimensions which arise due to the use of a UV cutoff as regulator 
and thus appear already at the perturbative level [26]. One possibility to handle these logarithmic 
divergences is to modify the integrand of the gluon loop such that the logarithmic divergences cancel 
[51, 54]. Another possibility is to employ a counter term, akin to the minimal subtraction scheme from 
perturbation theory [26, 27]. The divergences appear in the gluon two-point DSE in the following form: 



1 



Z(p 2 ) 



l + X(p 2 ) = l + 9 -Z'(p 2 ) + cg 



,lnA 2 



(5.3) 



where S(p 2 ) is logarithmically divergent and S'(p 2 ) is finite. Subtracting from this the self-energy at 
the momentum s times s 2 jp 2 yields 

11 



Z{p 2 ) Z{s 2 ) p 



P 



P 



P 



r£'(s 2 ), 



(5.4) 



i. e., the logarithmic divergence is gone. Since the logarithmic divergences only reside in the part pro- 
portional to p^Pv, its advantageous to use the longitudinally projected self-energy for the subtraction 
term since this interferes least with the physical part. The two methods for subtracting the logarithmic 
divergences are compared in Fig. 11. 



- 11 - 



G(p 2 ) Z(p 2 ) 




Figure 11. Ghost (left) and gluon (right) dressings obtained with a bare ghost-gluon vertex and the three- 
gluon vertex of Eq. (5.1) with a — 0.15. The red/solid lines result from subtracting the logarithmic divergences 
in the gluon loop kernel, while green/dashed lines stem from using the subtraction Eq. (5.4). 



5.1 Bare ghost-gluon vertex 

The use of a bare ghost-gluon vertex is standard in three and four dimensions, see, for example, 
[3, 4, 16, 17, 27, 30, 55]. The small deviations from the tree-level observed by lattice calculations 
[18, 56-58] do not affect the IR or UV behavior and influence only the mid-momentum regime. A 
semi-perturbative DSE analysis of the vertex [59] yielded qualitatively similar results as from lattice 
calculations. Therefore we start our investigations with a bare ghost-gluon vertex. The resulting 
propagator dressings are shown in Fig. 12. It turns out that the ghost dressing function does not 
approach one in the UV. However, since the dimension of the coupling constant is that of momentum 
we know from dimensional arguments that it should behave for large momenta as 

G(P 2 ) ^ ,+ \,2 (5-5) 
1 + cg z /p z 

with c a constant. Also lattice calculations [31] confirm that the ghost should be close to one in the 
UV. 

Based on this deviation of the ghost dressing function we conclude that the employed truncation 
is insufficient. The reason for this is a rather intricate relation in the ghost DSE. Its unsubtracted 
form has the structure 

1 + S(p 2 ). 



G(p 2 ) 

For large momenta, asymptotic freedom requires that S(p 2 ) vanishes, and thus the ghost dressing 
function approaches necessarily one. At low momenta, S(p 2 ) will behave in the scaling case as a + 
A(p 2 ) K + 0((p 2 ) 5 ) with some constant a and 6 > k. Only when a = —1, the scaling solution can be a 
solution to the equation. This requires that E must provide a large enough integrated strength. Due 
to the integral measure, the far infrared does not contribute to this, and due to asymptotic freedom 
neither does the ultraviolet. Thus, a = — 1 must be provided by the mid-momentum behavior. This 
can be generated by two ingredients. One is the mid-momentum enhancement of the gluon dressing, 
its maximum. The other is due to an enhancement of the ghost-gluon vertex. Here, we find that for a 
dressing function of the gluon, which is not exceeding the lattice results significantly, a = — 1 cannot 
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Figure 12. Propagator dressing functions using a bare ghost-gluon vertex and the three-gluon vertex ansatz 
given in Eq. (5.1) with a = 0.15,0.225. Left: Ghost dressing function. In the UV the values are 0.92 and 0.89. 
Right: Gluon dressing function. 



be achieved with a bare ghost-gluon vertex. This is different from three dimensions, where this is 
possible [27]. In four dimensions, this problem is alleviated, as the one in the ghost equation can be 
modified by a wave-function renormalization to any desired value [21], and the problem thus does not 
exist. 

5.2 Ghost-gluon vertex model 

Having identified the source of the problem of the UV behavior of the ghost dressing function we 
will now attempt to improve the results by using more elaborate ansatze for the three-point functions 
motivated by lattice results. In the following we will only use the transversely projected gluon DSE, 
so we only need to model the dressing function Df cc of the ghost gluon vertex which we take as 

D^oAr, P, q) = 1 + A2 ^ 2 \ 2 ^ 2 (fiR + fiM ^tfl ) ■ (5-6) 
A^ + p z + q z + r l \ A 4 + p 4 + g 4 J 

The new term proportional to Jim adds a bump in the mid-momentum regime as seen in lattice 
calculations [31]. This bump occurs also in three and four dimensions, see refs. [18, 56, 57], where it 
was also reproduced by a semi-perturbative DSE calculation [59]. The parameter sets the IR value 
of the dressing and A is a scale parameter. This form of the vertex leaves the value of n unchanged as 
discussed in Sec. 4.2, but not the IR value of the coupling, which scales with 1/(1 + fj^/A 2 ). 

For the three-gluon vertex we will continue to consider only one dressing given by Eq. (5.2), 
because adding dressings will not change anything qualitatively and one dressing is sufficient for our 
purposes. Fixing the parameters of the ghost-gluon vertex it is indeed possible to tune the three-gluon 
vertex parameters such that the ghost dressing function becomes one in the UV. A possible choice of 
parameter values is given in Table 1 and the resulting dressing functions are plotted in Fig. 13. 

A remarkable feature of these and the following results is that they cure the infrared singularities 
of perturbation theory. This is not so surprising when considering loops including gluons due to the 
infrared suppression of their propagator. For the ghost-loop in the gluon equation, however, this is 
naively unexpected, given that the ghost propagator is infrared more divergent as in perturbation 
theory. The explanation are intricate cancellations in the angular loop integral due to the non-trivial 
momentum dependencies of the dressing functions, emphasizing the delicate balancing in scaling type 
solutions. 
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Parameter 


Value 


A 


1 


fiR 


2.14189 


flM 


0.6 


hiR 


-71.5 


hiM 


9.88 



Table 1. Values employed for the parameters of the 
ghost-gluon and three-gluon vertex ansatze given in 
eqs. (5.6) and (5.2), respectively. 



Parameter 


Value 


A 


1 


hm 


-31.57 


hiM 


9.88 



Table 2. Values employed for the parameters of the 
three-gluon vertex ansatz given in Eq. (5.2) when in- 
cluding the ghost-gluon vertex dynamically. 
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Figure 13. Ghost (left) and gluon (right) dressings obtained with the models given in eqs. (5.6) and (5.2) 
using the values of Table 1 for the parameters. 



6 Including the ghost-gluon vertex dynamically 

Knowing how important the three-point functions are for the correct UV behavior of the ghost prop- 
agator we will now include the ghost-gluon vertex dynamically into our calculations. Working in 
two dimensions is thereby advantageous because we only have to calculate two-dimensional integrals, 
whereas in three and four dimensions the integrals are three-dimensional. Furthermore, in contrast 
to four dimensions, the UV behavior is trivial. This calculation is therefore also an exploratory study 
for future calculations in four dimensions which extend the currently employed truncation schemes 
beyond the propagators. 

We will again use the transversely projected gluon DSE so that only one dressing function of the 
ghost-gluon vertex is relevant. Three-point functions depend on three variables for which we choose the 
squares of the two ghost momenta and the angle ip between them. The vertex is calculated for a grid 
in these variables. For intermediate points we use linear interpolation. If any momentum is outside 
the grid we use the value of the dressing at its boundary. Considering the increased complexity of the 
ghost-gluon vertex DSE it was advantageous to derive the kernels with the program DoFun [44, 48]. 
For solving the DSEs the framework provided by CrasyDSE [60] was used. 

The ghost-gluon vertex has two distinct DSEs, which differ by the field that is attached to the bare 
vertex. We take the DSE where this is the gluon field. Although this DSE seems more complicated 
than the other one, since it has more terms, it turns out that in our truncation scheme it is simpler. 
The full DSE is shown in Fig. 14. The first step of our truncation consists in discarding all diagrams 
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Figure 14. The ghost-gluon vertex DSE. All internal propagators are dressed. Thick blobs denote dressed 
vertices. Wiggly lines are gluons, dashed ones ghosts. The employed truncation consists of the first, the second 
and the fourth diagrams. 



containing a bare four-gluon vertex, i. e., all two-loop diagrams. Of the remaining diagrams the ones 
with a bare ghost-gluon vertex dominate in the IR [28]. Furthermore, we keep only the diagrams of 
leading order in the UV, so we discard the third and fifth diagrams on the right-hand side. Although 
the latter is at leading order in the IR it does not represent one of the main contributions there, 
because when we insert the DSE of the irreducible quartic ghost vertex, we see that the IR leading 
contribution is a two-loop diagram and should give only a minor correction to the IR behavior of the 
ghost-gluon vertex. As attested by our calculations below already the fourth diagram yields only small 
corrections to the bare vertex. 

The resulting dressings of the propagators and the ghost-gluon vertex are shown in Figs. 15 and 
16, respectively. For comparison the plots of the propagator dressings also contain results presented 
in the previous sections. The differences originating in the mid-momentum regime can clearly be 
seen. Most notably the maximum of the gluon dressing function is driven to even larger momenta and 
becomes more shallow. In the ghost-gluon vertex dressing also a bump in the mid-momentum regime 
is seen. However, it can be uniquely traced back to the second diagram on the right-hand of its DSE 
in Fig. 14. We also tried the second version of the ghost-gluon vertex DSE, but we did not obtain a 
solution. The reason is that in this DSE the dressed instead of the bare three-gluon vertex appears 
and introduces an instability in the system of equations. This again illustrates the delicate balancing 
in the mid-momentum regime in two dimensions. 

Comparing the presented results with lattice data, we see that they do not agree as much as 
expected from three and four dimensions. Most likely these discrepancies could be compensated by 
the two-loop diagrams in the gluon equation and other neglected diagrams in the ghost-gluon vertex 
DSE. However, currently the inclusion of two-loop diagrams in DSEs has not been explored thoroughly 
enough to implement them here straightforwardly but see, for example, [61, 62]. The main reason 
that they have been neglected in the literature so far is that they yield only subleading contributions 
both in the IR and the UV for the Landau gauge in three and four dimensions [9, 28]. Furthermore, 
their increased complexity requires to extend currently employed methods [62]. Thus, though in two 
dimensions the qualitative properties of the solution of the DSEs is more simple and unambiguous, 



- 15 - 



G{ P Z ) 



Z(p z ) 





Z(p 2 )/p 2 



3.5 I *- 
I 



3.0 
2.5 
2.0 
1.5 

0.5 



S 

' \ 

/ s 



/ 
/ 

: r 
1* .- 



Figure 15. Ghost (left) and gluon (right) dressings obtained when including the ghost-gluon vertex dy- 
namically (red/solid lines) and using the three-gluon vertex ansatz from Eq. (5.2) with parameters given in 
Table 2 compared to the results when using a bare ghost-gluon vertex (green/dashed lines) or the model from 
Eq. (5.6) (blue/short-dashed line) with parameters given in Table 1. Bottom: Comparison of the propagators 
with lattice data (black points) from [26] . In the plot of the ghost propagators the continuum results are nearly 
indistinguishable. Scale set by matching the positions of the maxima in the gluon propagators. Note that this 
rescaling is responsible for the interchange of the blue/short-dashed line with the red/solid line. 



the quantitative features are not described as satisfactorily as in three and four dimensions. In four 
dimensions this difference is likely due to the possibility of rcnormalization, which permits to shift 
various effects to different momentum scales. In three dimensions, where one also finds quantitatively 
acceptable descriptions of the mid-momentum regime, it is probably due to the increased contribution 
from the momentum integral measure, which permits rather small differences at larger momenta to 
have already a significant impact. 

7 Conclusions 

Summarizing, we have provided the first full solution of the two-dimensional DSEs, including the 
equation for the propagators and the ghost-gluon vertex. The inclusion of the latter was necessary, as 
we found that the mixing of different momentum regimes in superrenormalizable theories invalidates 
the simple truncations used in the four-dimensional, renormalizable case. The reasons for this can 
be understood and lie in the intricate cancellations necessary for a scaling type solution, which we 
find, in accordance with [34], to be the only viable type of solution in two dimensions. That a similar 
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Figure 16. Dressing of the ghost-gluon vertex for various momentum configurations obtained from the 
coupled system of propagators and the vertex itself. For the three-gluon vertex the ansatz of Eq. (5.2) with 
the parameters of Table 2 was used. Top: Fixed angle as indicated at the top of each plot. Middle: Fixed 
momentum q 2 as indicated at the top of each plot. Note the different scales on the z-axes. Bottom: On 
the left the symmetric configuration is shown. Since this configuration cannot be realized on the lattice, no 
comparison with lattice data is possible. The bump stems from the non-Abelian diagram. On the right the 
gluon momentum and one ghost momentum are orthogonal. The red/solid line is the result from the DSE 
calculation and green squares are for = 10/L = 21/m" 1 and blue diamonds for /3 = 22.5/i = 12/m -1 . 
Lattice data from Ref. [31]. 
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situation did not arise in this severeness in earlier studies in three dimensions [27] must be attributed 
with hindsight to quantitative effects. 

As a consequence, the results here have not yet reached a maturity as in higher dimensions when 
it comes to quantitatively reproducing lattice results. This will require much more sophisticated trun- 
cation schemes, which will likely require the inclusion of two-loop terms. In DSEs, this is a formidable 
endeavor, and thus a renormalization group approach, with its intrinsic one-loop structure, may be 
more adequate. Nonetheless, we reproduced all qualitative features and provided an understanding of 
the underlying mechanisms. 
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A Kernels 

The kernel K G of the ghost DSE given in Eq. (2.7) is 



(x 2 + (y - z) 2 - 2 x(y + zj) 
Axy 2 z 



K G{p,q) = - jz^jz (A- 1 ) 



with x =p 2 , y = q 2 and z = (p + q) 2 . The two kernels K 9 z h and K 9 J of the gluon DSE Eq. (2.8) read 
K 9 z h (p,q) 



,,,/,, n n] = {x 2 + {y - z) 2 - 2x{y + z)) 
4x 2 yz 



t x i ~8xyz(y + z)+x 2 (-2y 2 -8yz~2z 2 ) + (y-z) 2 (y 2 + 2yz + z 2 ) 

K Z (P> 1) = " o 9 9 9 • ( A - 2 ) 



8x 2 y 2 z 2 

In order to get rid of the logarithmic divergences in the gluon DSE without using counter terms, the 

rgl. 



following expression is added to the kernel of the gluon loop K 9 z l - 



K z ,8Uh M = ^ y - (A.3) 

The kernels for the ghost-gluon vertex are too lengthy to be reproduced here. They were generated 
automatically using the programs DoFun [44, 48] and CrasyDSE [60] using the ansatze for the vertices 
as described in the main text. 



B Three-gluon vertex 

For solving the coupled system of propagators and ghost-gluon vertex in Sec. 6 we employed a lattice 
inspired ansatz for the three-gluon vertex. Here we calculate the three-gluon vertex using the results 
obtained there for illustration and comparison. A combined calculation of propagators and both 
three-point functions did not yield a stable iteration which is again due to the effect truncations have 
on the mid-momentum regime. Nevertheless it is interesting to see that the lattice results can be 
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Figure 17. The three-gluon vertex DSE. All internal propagators are dressed. Thick blobs denote dressed 
vertices. Wiggly lines are gluons, dashed ones ghosts. 



reproduced at least qualitatively. Whereas the correct IR behavior should emerge automatically, the 
most important question is if the zero crossing as observed on the lattice is reproduced [31]. 

The DSE for the three-gluon vertex is depicted in Fig. 17. In analogy to the truncation for the 
ghost-gluon vertex we discarded all diagrams with two-loops or non-primitive vertices. This provides 
the correct IR behavior and includes also the leading corrections to the tree-level behavior in the UV. 
The remaining loops are the so-called ghost and gluon triangles represented by the second and fourth 
diagrams, respectively, on the right-hand side in Fig. 17 and the swordfish diagrams represented by 
the fifth, sixth and seventh diagrams. For the dressed four-gluon vertex we use the bare one. This 
does not affect the IR and UV behavior but only the mid-momentum regime. Another approximation 
is that we use the result for the dressing of the three-gluon vertex as given in Eq. (B.l) below on the 
right-hand side of the three-gluon vertex DSE, i. e., T A = D^ ro y 

In lattice calculations the calculated quantity is the contraction of the transversely projected 
three-gluon vertex with the bare three-gluon vertex [31, 57] normalized to the tree-level: 



This quantity was calculated using a standard fixed-point iteration. The results for some selected 
momentum configurations are shown in Fig. 18. The qualitative features expected from lattice cal- 
culations, especially the zero crossing, are all there. At the symmetric point we compare the results 
when taking into account only the ghost triangle, both triangles or all five diagrams of the trunca- 
tion described above. For the so-called orthogonal momentum configuration a comparison between 
continuum and lattice results is shown in Fig. 19, where again also the influence of the gluon triangle 
and the swordfish diagrams is shown. As expected the gluon triangle only gives a contribution in the 
mid- momentum regime and creates a small bump there which is absent when only the ghost triangle is 
taken into account. Interestingly, the swordfish diagrams reverse that effect and the bump goes away 
again. The zero crossing is also almost unaffected by the inclusion of further diagrams. Note that 
the small influence of gluons in the three-gluon vertex DSE is a non-trivial result insofar as for the 
ghost-triangle-only truncation the three-gluon vertex does not appear on the right-hand side and the 
iteration consists of only one step. Only the inclusion of the gluonic diagrams leads to the feedback of 
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Figure 18. Dressing of the three-gluon vertex, see Eq. (B.l). Only the beginning of the IR divergence is 
shown by cutting all data below —10. Top left: Angle fixed at arccos(— 0.41625). Top right: One momentum is 
fixed at V0. 00007786 g. Bottom right: Symmetric configuration. The solid/red line is with the ghost triangle 
only, the dashed/green one with both triangles and the dotted/blue line with all five diagrams. 



the three-gluon vertex onto itself. As expected, the details of the mid-momentum regime depend on 
the employed truncation, whereas the IR and UV regime are unaffected. 
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